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ABSTRACT 

The paper describes a new technique for convert- 
ing a constrained optimization problem to an 
unconstrained one, and a new method for multi - 
objective optimization based on that technique. 
The technique transforms the objective functions 
Into goal constraints. The goal constraints are 
appended to the set of behavior constraints and 
the envelope of all functions In the set Is 
searched for an unconstrained minimum. The 
technique may be categorized as a SUMT algor- 
ithm. In multi -objective applications, the 
approach has the advantage of locating a com- 
promise minimum without the need to optimize for 
each Individual objective function separately. 
The constrained to unconstrained conversion Is 
described, followed by a description of the 
multiobjective problem. Two example problems 
are presented to demonstrate the robustness of 
the method. 


max maximum value 

min minimum value 

INTRODUCTION 

The paper has two purposes: to bring to the 

optimization practitioners* attention a new 
technique for converting a constrained optimiz- 
ation problem to an unconstrained one, and to 
present a new method for multiobjective optimiz- 
ation based on that technique. 

The conversion technique 1 may be categorized as 
a "Sequential Unconstrained Minimization Techni- 
que" (SUMT) class method 2 but It does not 
require the use of a draw-down factor, unlike 
the classical procedure. Also, the unconstrain- 
ed function It uses to represent the constrained 
problem at hand Is defined over both the feas- 
ible and Infeasible domains, similar to an 
extended penalty function . 


NOMENCLATURE 

A cross sectional area 

E 4 Youngs Modulus of member 1 

F|< k-th objective function 

Fe 1 ron 

F global criterion or compromise objec- 

tive function 

F* reduced objective function or goal 

constraint 

fj set of functions 

problem constraints 
J NCON + NOBJ 

K-S Krelsselmeler-Stelnhauser function 

NCON number of constraints 

NOBJ number of objective functions 

T1 titanium 

X vector of design variables 

p K-$ coefficient 

Superscripts 


Difficulty in defining a single objective func- 
tion In many engineering design problems Is a 
motivation for continuing Interest In develop- 
ment of techniques for multi objective optimiza- 
tion applications. Many of the mul tl object 1 ve 
optimization methods require either a conversion 
to a single objective function by means of a 
composite function with judgmental "weight fac- 
tors", or separate optimizations for each objec- 
tive followed by an additional "global" optimiz- 
ation to arrive at a suitable compromise . The 
technique Introduced In this paper, for the con- 
stralned-to-unconstralned optimization problem 
conversion. Is also shown to have an Intrinsic 
applicability to multiobjective optimization. 

Its primary benefit In that application Is the 
elimination of the potentially expensive sepa- 
rate optimizations for each objective. 

CONSTRAINED-TO-UNCON$TRAINED CONVERSION 


cm constrained minimum 

0 Initial design point 

s scaled value 

Subscripts 

1 1 th element 

j jth element 

k kth element 
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The conversion technique replaces the constraint 
boundary surfaces and the objective function 
surface In n-dlmenslonal design space with a 
single envelope surface constructed uslno the 
Krelsselmeler-Stelnhauser (K-S) function* first 
Introduced for optimization of control systems. 
That function has subsequently been established 
In structural optimization as a means to replace 
many constraints with a single cumulative con- 
straint 6 . The function Is a differentiable 
envelope of a set of functions f j (X) and It 
has this, form: 


K-S 


f max + 
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and a property such that 
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A Log(J) 
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( 2 ) 


where p control* the distance of the K-S enve- 
lope surface froai the fjpx surface. The K-$ 
function may be regarded** analogous to the MAX 
function available In many high-order program- 
ming languages but» unlike the MAX function. It 
Is differentiable (value- and slope-continuous); 
therefore, It may be called upon by a gradient- 
guided optimization algorithm to search for a 
minimum of the envelope of a set of functions. 

The constrained optimization problem to be 
solved with the aid of the K-S function Is, In 
conventional formulation: 

minimize F k (X), k * 1 to NOBJ; (3) 


such that 

gi (X) <0, 1 » 1 to NOON; 

where constraint functions g^ are written In 
terms of the computable functions, termed DEMAND 
(X) and CAPACITY (X), that provide the measures, 
respectively, of what the design Is asked to 
carry versus what It can sustain: 

gi (X) * DEMAND ( X ) /CAPAC I TY ( X ) -1 (4) 

For Introductory purposes, F k (X) in equation 3 
Is a single objective (l.e., NOBJ * 1); exten- 
sion to many objectives will follow later. To 
formulate the K-S function as an envelope of the 
objective functions and constraints, one has to 
normalize the objective function In order to 
make It comparable to the normalized constraint 
functions. The normalized objective and con- 
straint functions form a set of functions whose 
envelope Is approximated by the K-S function. 

An unconstrained minimum (except for the usual 
side constraints) of the K-S envelope may be 
found by any suitable search algorithm. 

The procedure formulated 1 Is Illustrated in 
figures 1 and 2. For graphic simplicity, one 
design variable 1$ shown. To make this paper 
self-contained the procedure Is restated here In 
descriptive terms keyed to figures 1 and 2. 

A single design variable x Is measured on the 
horizontal axis. The objective function F and 
contraint functions gi and g 2 are represented on 
the vertical axis. The Initial design point Is 
at x « x° where the constraints are violated and 
the objective has the value F*. By inspection, 
the constrained minimum lies at x cm . The 
requirement Is to locate that minimum starting 
from x°. 

Referring to figure 1, the objective function In 
Its original form before normalization Is 
labeled F. F is normalized by dividing F by 
F°. The scaled F s * F/jF°| then shifted to 
Intersect the abscissa by subtracting unity and 


further offset by subtracting g^x. The 
shifting and offsetting Is expressed as follows: 

*( F s - 1) - «mx W 

which moves the objective function to F* In 
figure 1. The normalized, shifted, and offset 
objective function, F will be referred to as 
a reduced objective function. The reduced 
objective function Is Included with the con- 
straint functions to form a set of functions 
fj(x) whose envelope Is approximated by the 
K-S function shown by the dashed line. 

An unconstrained minimum of that K-S function Is 
found at x 1 by means of any search method suit- 
able for unconstrained optimization. Locating 
that minimum completes one cycle of the proce- 
dure. The procedure cycle count should not be 
confused with the count of the Iterations carri- 
ed out by the unconstrained minimum search algo- 
rithm; many of the latter Iterations may be 
needed In one of the former. 

Referring to figure 2, the next cycle starts 
with x l and F . Equation 5 Is used to compute a 
new F* using F l . This formulation takes Into 
account that g^x may be a negative value. 

The K-S function Is fitted to the set (gi ,$ ,F*) 
using equation 1 and Its minimum Is found at x . 
This completes cycle 2 of the optimization pro- 
cedure. Successive cycles are continued until 
convergence. 

It Is apparent, from the above two cycles, that 
the unconstrained minima progress from the 
Initial x® to x , to x , approaching the con- 
strained minimum at x cm . One may also observe 
that In contrast to cycle 1, cycle 2 starts from 
a feasible design. This illustrates the capa- 
bility of proceeding either from an Infeasible 
or a feasible Initial design point. Finally, It 
should be noted that the process of shifting and 
offsetting changes the position of the normaliz- 
ed objective but perserves Its slope for each 
cycle. 

At the Initial location x°, the K-S function 
reflects almost exclusively the geometrical pro- 
perties of the constraint boundary gi . In con- 
trast, the objective function dominates the K-S 
envelope function In the search for a minimum In 
the case of a feasible design. Consequently , If 
a design point In the midst of the 1-th cycle Is 
Infeasible, the search direction toward a small- 
er K-S envelope will point toward smaller values 
of the dominant constraint, thus reducing the 
amount of constraint violation at the possible 
expense of Increasing the objective function. 

On the other hand, If that design Is feasible, 
the search direction toward reduced K-S values 
will be equivalent to moving toward lesser 
values of the objective at the possible expense 
of increasing the values of the satisfied domi- 
nant constraints. 

The process converges to the state shown in 
figure 3 when It Is no longer possible to reduce 
the objective function without violating the 
constraints. At the constrained minimum point 
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the formula of equation 5 produces no offset 
because g^x * 0 . The constrained mini mum 1 $ 
properly located at the boundary surface of the 
maximal constraint. Generalized to n-dlmenslon- 
al space with m constraints, the unconstrained 
minimum of the K-S envelope constructed as above 
approximates the location of a constrained mini- 
mum defined by • full vertex of the design 
space, or by « point of tangency between the 
objective function and the dominant constraint, 
or a mixture of these two extreme conditions. 


MULTIOBJECTIVE OPTIMIZATION 


The same procedure described above, applied In 
this case, will establish a sequence of K-$ 
envelope minima with each K-S envelope function 
containing all of the reduced objective func- 
tions. Since all normalized objective functions 
Intersect at the same point located at x°, their 
values In the Interval x > x° Into which the 
search will be progressing will rank according 
to the magnitudes of the corresponding slopes 
(for a convex problem). The sequence leads to a 
constrained minimum point, where the configura- 
tion of constraint boundaries and reduced objec- 
tive functions are Illustrated In figure 4. The 
configuration Is Identical to that shown In fig- 
ure 3 except for the presence of several objec- 
tive functions, all of which Intersect at the 
constrained minimum point after scaling and off- 
setting. The constrained multlobjectl ve minimum 
point, Illustrated In figure 4, has the property 
fthat one can not depart from It without either 
violating the constraint (s) or Increasing at 
least one of the objective functions - the clas- 
sical definition of a pareto-optlmum. It Is 
pointed out In the appendix how the method 
relates to the goal programming class of algor- 
ithms. 


NUMERICAL EXAMPLES 


The method was tested using a three-bar truss 
design problem, figure 5, which Is a classical 
test case . The truss Is considered symmetric, 
thus Ai ■ As and £\ * E 3 . This problem consi- 
ders two distinct load cases Pi and P 2 . The 
problem Is generalized by Including the angle a 
with the cross-sectional areas Ai , A 2 as a third 
design variable. In addition, two materials are 
used. Ai and As are made of the same material 
and Aa Is constructed of a different material. 
The goal of the optimization is to minimize 
combinations of weight, cost and support area 
width 0 In figure 5. Steel, with a greater 
Youngs Modulus and Inexpensive, and titanium, 
with a greater stress allowable, a lower Youngs 
Modulus, and much more expensive were selected 
as the two materials to provide for a meaningful 
trade-off between the objectives of weight and 
cost. Introduction of the angle o brings In the 
dimension D, a function of <*, as a consideration 
In design; reduction of D Increases the forces 
In the members and also reduces the length of 
the outer members, such that a complex coupling 
of the geometry, weight, cost, and strength Is 
created. The truss analysis Is given In the 


Appendix and the material properties. Including 
costs, are stated In table 1 . 

Verification of the method with one material 
common to all the members, the objective of min- 
imum weight, the angle « fixed at 45 degrees, 
and the two cross-sectional areas Ai and Aa as 
the only two design variables, yielded results 
which agreed with the classical test case . 
Subsequently, a number of optimization experi- 
ments were carried out with various combinations 
of material configurations, choice of design 
variables, and selection of objectives. In all 
cases the search for the minimum of the K-S 
function envelope was carried out g by the 
Oavldon-Fletcher-Powell algorithm . 

Table 2A displays a list of cases, showing for 
each case the objective functlon(s), material (s) 
used for the truss members and design vari- 
able^). For example, case 9 In table 2A uses 
weight and the dimension D as objective func- 
tions; steel for the outboard members; titanium 
for the center member; and Ai , A 2 , and a as 
design variables. Initial conditions of Ai , Aa , 
and a for all cases of table 2A are given In 
parenthesis, table 2B. 

The results In table 2B, corresponding to the 
cases of table 2A, show consistently the 
method’s ability to generate optima for single 
and several objectives for the example problem. 
It Is Instructive to observe that when an opti- 
mization Is executed with a single objective, 
for example, the weight, and then repeated with 
additional objectives, the latter are signifi- 
cantly reduced at the price of a relatively 
small Increase of the former - see for Instance 
the cases 10 , 12 , and 16 for (WFa/Tj) 
and cases 9, 11, and 13 for (F e /Tj/F e ). 

These cases are also Illustrated graphically In 
figures 6 and 7. 


Cases 15 and 18 In table 2B Include results 
obtained by the global criterion formulation . 
That method solves a multlobjectl ve optimization 
problem by first executing separate optimiza- 
tions for each of the objectives F^ to obtain 
a set of feasible minimum solutions Ft.. 

Next, a search Is carried out for a con- 
strained minimum of a compromise objective F(X) 
where 


F(X) 


[M*) h min y 

NOBJ 

k*l ^ k m 1 n 


( 6 ) 


such that 

91 (X) < 0, 1 ■ 1 to NCON; 

Comparison with cases 19 and 20 respectively 
shows a good match between the method reported 
herein and the global method, figure 7. 

Additional test cases, not Included In the 
above, were carried oat starting the optimiza- 
tion procedure from a variety of Initially feas- 
ible and Infeasible points. Convergence to the 


ORIGINAL PAGE IS 

3 QE EOOR QUALITY 



same optimal design was observed, regardless of 
the Initial point, except for cases which 
clearly Indicated existence of local minima. 

One should note that by adjusting the 
coefficient p In equation 1 toward smaller 
values, one can keep the K-S function from 
•dipping" Into Inferior local minima, thus 
maintaining the search advance toward the global 
minimum. A method for controlling p has been 
suggested . Typically p may be held constant 
after a proper value Is determined. It should 
be noted that the coefficient p affords the user 
a degree of control over the procedure. On the 
other hand, problem dependent sensitivity of the 
procedure to that coefficient may be a drawback 
which forces experimentation to find a range of 
P values for the case at hand. 

The second example Is an optimization of an 
electric power transmission line shown In figure 
8. The nomenclature for this problem Is listed 
In the appendix. The line Is assumed to extend 
over a fixed distance D, over a flat terrain. 

It comprises equidistant towers, separated by a 
distance L, assumed to be thin-walled cylindri- 
cal columns of radius R, with a constant wall 
thickness t. All towers are geometrically iden- 
tical and made of steel. Differences between 
the weight and cost of the real towers and their 
Idealization as cylindrical columns are repre- 
sented by shape factors. The towers support 
three parallel electric cables made of an alumi- 
num alloy. The cable weight and the weight of 
Ice accumulated on the cables are the combined 
ioads q, which put the cables In tension and the 
towers In compression. The cables are assumed 
to be Infinitely rigid In tension so that they 
form a catenary curve between the towers. The 
catenary sag Is assumed small relative to the 
cable length so that a parabolic approximation 
to the catenary equation Is used. Constraints 
are placed on the stresses In the towers due to 
yielding, cylinder wall buckling, and column 
buckling. Coluim bending Is also considered for 
the case of failure of all cables at a point In 
the line. Tensile stress constraints are Impos- 
ed for the electric cables, and the cable-to- 
ground clearance Is constrained. Side con- 
straints are Imposed on the cable cross-section- 
al areas A, the tower geometry, and the distance 
between towers L. Another side constraint Is 
prescribed for a minimal cable cross sectional 
area A as required for electric power transmis- 
sion. Problem design variables are the distance 
between towers L, tower height H, column radius 
R, wall thickness t, cable cross-sectional area 
A, and cable tension, P. The problem objectives 
are the total weight of steel and aluminum W, 
the total cost, and the tower height. The cost 
objective entails the cost of: steel, aluminum, 

and tower foundations. The height objective Is 
Included to reflect environmental specifica- 
tions, a typical unquantlflable consideration 
engineers are being confronted with Increasingly 
often. 

This case Is rich In complex Interaction among 
the design variables. For example, the required 
cable-to-ground clearance may be attained by at 
least three different means. Increasing the 
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tower height, cable tension, or decreasing the 
tower-to-tower distance. However, each of these 
means has a different Influence on each of the 
objective functions and constraints. 

Details of the example analysis are provided In 
the Appendix. Table 3A displays a list of cases 
for this example. Optimization results for com- 
binations of the objective functions are given 
In table 3B. Initial conditions of design vari- 
ables for all cases are given In parenthesis, 
table 3B. Figure 9 displays results of several 
cases of discussion graphically. The results, 
again, show that the method has the ability to 
locate compromise designs satisfying all the 
constraints. The optimization for cost only 
(case 2) reduces the cost by 28% of the cost of 
the minimum weight design, case 1, while Increa- 
sing the weight by only 9X. In contrast, the 
reduction of the tower height objective is much 
more expensive In terms of the weight and cost 
(because the cable tension tends to infinity as 
the tower height approaches the required cable 
ground clearance). Thus, the inclusion of the 
tower height as another objective, along with 
weight and cost, results in the cost, case 6, 
nearly 12 times greater and weight nearly 3 
times greater than the weight and cost of the 
minimum cost only design, case 2. In the latter 
case, the towers are widely separated and made 
tall which reduces the cost of the tower founda- 
tions. On the other hand, the cost and weight 
of the compromise design, case 6, are 7% and BOX 
lower than those of the design for minimum 
height alone, case 3. 

CONCLUDING REMARKS 

A new technique 1 for converting a constrained 
minima problem to an unconstrained one was 
demonstrated to be a useful tool In single 
objective and multi objective applications. The 
technique transforms the objective functions 
Into goal constraints, the goal value for each 
objective Is an adjustable quantity. The objec- 
tive goal constraints are then appended to the 
set of behavior constraints and the envelope 
(cumulative constraint) to all the functions In 
the set Is constructed using the Krelsselmeler- 
Stelnhauser function, whose minimum Is searched 
for by any unconstrained minimization algorithm 
(the Davidon-Fletcher-Powell method 8 was used In 
the reported study). Search toward the minimum 
of the envelope function advances the design 
toward the compromise constrained minimum. That 
minimum Is reached In an Iterative procedure, 
which updates a set of behavior and objective 
goal constraints and their envelope function at 
the outset of each Iterative cycle. By repre- 
senting the objective functlon(s) as objective 
goal constraints, the method Is related to the 
goal programming approach , and the constrained 
minimum It attains conforms to the classical 
pareto -optimum definition. 

The method typically converged after 8 to 50 
cycles, depending on the mix of design vari- 
ables, parameters andobjectlve functions. 

The technique was demonstrated on variable 
geometry and cross-section trusses built of 
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different mrterlals contributing a wide range of 
Mechanical properties and cost to the design 
objectives. Objectives also Included the amount 
of space occupied by the truss on Its support 
surface. An example of an electric power 
transmission line was also optimized for a 
compromise of objectives . The objectives were 
material voTwe* «e«t* and support tower height. 

The method results compared well with those 
obtained by a goal programming algorithm and the 
method performance was satisfactory In all of 
the single and multi objective test problems. In 
contrast to other multi objective optimization 
procedures, the method showed an ability to 
locate compromise optimum designs without the 
expense of having to optimize individual 
objectives. 
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APPENDIX 

Step-by-Step Optimization Procedure 

1. Set Iteration counter K-0; Initialize x • 
x*; 

2. Execute analysis of the problem to 
obtain Ffc and g^'s; 

3. Define a reduced objective function per 
equation 5; 

4. Define an envelope function K-S(F*,gj) per 
equation 1; 

5. Find coordinates X . of the 
unconstrained (except for side constraints) 
minimum of the K-S function by any 
sultableunconstralned optimization 
algorithm. 

6. Reset K * K + 1, and reset x* * x ^ n ; 

7. Repeat from 12 until convergence criteria Is 
satisfied. 

A computational cost saving option In step 5 Is 
to execute only a limited number of steps toward 
the K-S minimum Instead of a full fledged func- 
tion minimization. Under that option, one may 
begin with a single step and progressively 
Increase the number of steps as the procedure 
advances. In the case of many objective func- 
tions step 3 Is carried out for all of thej, and 
the entire set of the reduced objectives F 
are included as arguments In the K-S function In 
step 4. 

Analysis of the Truss Example 

Referring to figure 5, the truss Is analyzed as 
a system with two elastic degrees of freedom. 
Utilizing the stiffness method the following 
load/deflection equations are obtained: 
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Kn 

K 12 

K 21 

*22 


ui 

U2 


(Ril: {Ra) ; 


Due to symmetry K 12 9 Kzi . 
Where: 


Kn * cos 2 o(k3 ♦ ki ); 

Ki 2 a -cosa k3 slna + cosa ki slna; 


9. Barthelemy, J-F«t M. Riley, M.F.: "An K 22 * s1n 2 o(k3 + ki ) + k 2 ; 

Improved Multilevel Optimization Approach 
for the Design of Complex Engineering 
Systems", AIAA paper 86-0950-CP, presented and: 

at the AIAA/ASME/ASCE/AHS 27th Structures, 

Structural Dynamics and Materials DET * K 11 K 22 - K: ; 

Conference, San Antonio, TX, May 1986. . 2 
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The member stiffness are given as, 

kj * EjAj/a j 9 J * 1,.*.3 

where: 

t\ • h/slnm 
i2 * h 
is - ij 


For load case Pi : 


Ri 


-Pi cose 
Pi sine 


For load case P2 : 

I P2 cos8 

; 

P2 slnB 


The solution of equation A1 for load case Pi 
yields: 

ui * DET“*((-P cos® ) k2 2 - (Pi s1na)ki2); 

u 2 « DET' 1 ((Pi slna)kn - (-P cosajku); 

and for load case P2 yields: 

1 mi « DET ^ 1 ( (P2 cos®)k22 - (P2 sin® )ki 2 ) ; 

U2 * 0ET“ 1 ((P2 $1n®)kn - (P2 cosa)ki2); A3 


From the displacements A3 and Young’s moduli 
Ea for each member the stresses are recovered 
a* follows: . 


H - tower height; 

0 - total distance covered by the transmission 

line (60.0 ml.); 

1 - true length of the cable between towers; 

A - cable cross-section area (minimum 0.20 

In/); 

P - cable tension; 

f - sag of the catenary between towers; 

E - Young’s modulus for steel (towers) (3.0 x 
1 <f psl.); 

ri - specific weight of aluminum (0.100 

lb/in. 3 ); 3 

Y2 - specific weight of steel (0.284 lb. /in. ); 
Y3 - specific weight of Ice (0.033 lb/1 n . ); 
mi - shape factor for tower weight (0.75); 
to - shape factor for tower cost (1.0); 
n - number of parallel cables; (3); 
so - safety factor for cables (1.5); 
si - safety factor for combined bending- 

compression stress In the tower (1.25); 

52 - safety factor for cable stress (1.50); 

53 - safety factor for compressive stress In 

tower (2.00); 

s n - buckling factor (0.60); 
q - load on the cable Including cable weight and 
ice accumulation; 

ci - cost of the cable material per unit weight 
(0.40 J/lb.); 

C2 - cost of the tower material per unit weight 
(0.09 $/lb. ) ; 

C3 - cost of the foundation per tower 
($50,000.00); 

k voi - ratio of the volume of Ice to the 
volume of the cable (3); 
o 1 - allowable stress for the cable (44 ksl); 

®2 - allowable stress for the tower (36 ksl); 
h - cable-to-ground clearance (H-f); 
h P - required cable-to-ground clearance (32.8 
ft.); 

The strength of materials equations for the 
problem are: 


®i * ((ui cosa + U2 slna) / ii )Ei ; 


Weight on the unit length of the cable Including 
Ice 


o 2 « (U2 /i2 )&2 *» 


®3 ■ ((-UI cosa + U 2 slna) / ij )Ej ; A4 


Analysis of the Electric Power Transmission Line 
Example 

Referlng to figure 8, the material used for the 
cables Is 2024 aluminum alloy. The material 
used for the towers Is AISI carbon steel. The 
following Is the nomenclature used In the analy- 
sis. The actual data are noted In parentheses. 
Including minimum (lower bounds) values where 
appropriate. 

R - tower mean cylindrical cross-section radius 
(minimum 5.0 In.); 

L - tower -to -tower distance (minimum 300.0 
In.); 

t - tower wall thickness (minimum 0.125 In.); 
bi - ratio of R/t (minimum 50); 
b2 - ratio of l/H (minimum 4); 


q * Y2A(1.0 + k (Y 3 /Y 2 )) A5 

vOl 

Sag of the cable and Its true length (assuming 
Inextenslonal cable) 

f = q L 2 / (8 P) A6 

1 * 1(1.0 + 2.67(f Z /L)) A7 

Stress In the cable 

0 = $2 P/A A8 

c 


Stress In the tower due to combined compression 
and bending; the worst case of the latter occurs 
when the cables break on one side of the tower 

° T * (nPH/»R 2 t) ♦ (nP(4f/L)) A9 
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Critical stress for cylindrical wall buckling 
(conservatively using a formula for uniform 
compression although the tower Is in the state 
of combined compression and bending) 

a * t| 8 f n f/L A10 

comp 

Critical force of the tower buckling In the 
column mode 

F - « 3 ER 3 t/(4 H 2 ) All 

crlt 

The total material weight of the towers and 
cables (neglecting the last tower) 

W t « (Alnyi + mi 2 »RtHY 2 )(D/L) A12 

The total cost of the transmission line and 
towers Including the foundation costs 

Cost - ((ciAlTin(l.O + 2.67(f 2 /L) )) + 


(ns2wc2RtHt2) ♦ C3)(D/L) A13 

The objectives of the problem are contributed by 
eq. A12, A13, and the tower height H. The con- 
straints are: 

cable to ground clearance 

* > h > h r A14 

strength of the tower 

< fr 2 /si A15 

buckling of the tower In the cylindrical wall 
and column modes* respectively 


°comp < °cy! 2 *^/s2 

A16 

where: ® cy1 - &.E(t/R)/s 2 

°comp < f CRIT/*3 

AI7 

strength of the cable 

®c £«»/*<> 

A18 

Ratio of R/t 

R/t > bi 

A19 

Ratio of L/H 


l/H > t» 

A20 

The design variables are A, R, t. 

H* L and P 


The Method's Relationship to Goal Programming 

It can be shown that the method described 
In the report relates to the methods of the goal 
programming category. Indeed* by Including 
coefficients A and n In equation 5 


F k A21 

F\ * ( — " n k) ' %ax 

|p%| 6 k 

where: F% are previously computed feasable 
solutions, It becomes possible to assign prior- 
ities to the objective functions by controlling 
the shifted distance and slopes (l.e., the 
relative magnitudes of the normalized objectives 
away from their Intersection point). One may 
observe that at the multi objective constrained 
minimum point we have 


- * 0 A22 


for all the objective functions Fj<. By defi- 
nition of F|< , the expression In parentheses 
may be regarded as a constraint Imposed on the 
objective F%. If the value F°k m1n Is 

known, one could replace F^^ for |^°k| 

and 1 for In equation A22 and, thus, 
formulate the optimization problem as follows. 

Find X such that 

Fk/F|c 4 - 1 < 0; k » 1*...N0BJ 

* *min — 


gi i 0, 1 * 1....NC0N; 


This formulation states a solution to a set of 
Inequalities* expressing the Intent to modify 
the design to a state In which all the behavior- 
al constraints are satisfied and all objectives 
are maintained below or at their target 
linin’ Thus procedure may also be 
categorized as a goal programming approach. 
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Table 1. Material Properties and Coats 



HAT El I At | 

STta (Ft) 

TITANIUM (T1) 

pacific 

i.nt 

0.110 

Cast (Vibe) 

0.41 

35.00 

A11a«01t itrtu-t«nt<w (til) 

n 

110 

itrm^iNprtulw (fcll) 

37 

83.5 

rovsf'i nadwivs (mD 

30 i 10 1 

15.5 i 10* 


# 


Table 2A. Cases for Three-Bar Truss Problem Table 2B. Results for Three-Bar Truss Problem 


1 CASE 

WIGHT (lb) 

COST (S) 

Ai (In*) 

At (In*) 

• (dtj) 

Initial Condition* of Dtslgn Variables 

(1.0) 

O-o) 

(45.0) 

1 

4.18 

1.71 

0.441 

0.3341 

45.00 

3 

4.30 

17.56 

0.450 

0.403 

45.00 

3 

0.805 

14.05 

0.133 

0.088 

45.00 

4 

4.54 

3.34 

0.563 

0.013 

45.00 

S 

4.43 

3.45 

p.548 

0.011 

45.00 

4 

1.09 

13.31 

.108 

.311 

45.00 

7 

0.936 

13.44 

.106 

.301 

45.00 

8 

4.091 

3.50 

i«003 

0.039 

75.44 

9 

4.333 

5.33 

1.074 

0.067 

76.89 

10 

3.7« 

63.30 

0.779 

0.453 

•4.03 

u 

4.881 

8.640 

0.767 

0.169 

69.67 

1? 

1.153 

14.69 

0.157 

0.304 

80.15 

13 

4.117 

10.30 

0.536 

0.3161 

53.35 

14 

4.313 

3.19 

0.433 

0.013 

58.37 

15 

4.148 

3.13 

0.403 

0.011 

55.34 

14 

0.844 

17.89 

0.199 

0.054 

63.65 

17 

1.093 

13.30 

0.105 

0.318 

44.43 

18 

1.143 

13.30 

0.103 

0.336 

43.79 

19 

4.004 

3.05 

0.593 

0.010 

54.75 

30 

0.984 

13.43 

0.103 

0.176 

43.11 


CASE 

08JECTIVC 

FUNCTION(S) 

WTE4IAUS) 

DESIGN 

VAA/4fl£(S) 

1 

Ualfht 

Fa/Ft/Ft 

A 

3 


Fa/T1/Fa 

A 

3 

tfclgfct 

T1/F*/T1 

A 

4 

Cost 

Fa/T1/Fa 

A 

5 

*tV*t/Cost 

Fa/T1/Fa 

A 

6 

Cost 

TI/Ff/TI 

A 

7 

Halgftt/Cast 

TI/Ft/TI 

A 

8 

Ualfht/D 

Fa/Fa/Ft 

Va 

9 

HH0ht/D 

Ft/T1/Fa 

v. 

10 

Halflit/0 

TI/Fa/TI 

A/« 

11 

W /Cost/D 

Ft/T1/Ft 

A/a 

13 

m/cost/D 

TI/Fa/TI 

A/a 

13 

Htlfht 

Ft/Tl/Fa 

Va 

14 

Cost 

F*/T1/Ft 

Va 

IS 

Goal 

Fa/T1/Fa 

Va 

16 

Malfftt 

TI/Fa/TI 

Va 

17 

Cost 

TI/Fa/TI 

A/a 

18 

OMl 

TI/Fa/TI 

A/a 

19 

mifAt/Coit 

Fa/TI/Fa 

Va 

30 

Hal 98t /Cost 

TI/Fa/TI 

Va 


Table 3 A . Cases for Transmission 
Tower and Cables 


Table 3B. Results for Transmission Tower and Cables 


CASE 

Halaht (iba.) 
(HO*) 

Cost (!) 
(110*) 

9 (In.) 

H (IB.) 

t (I*.) 

« d».») 

f 08.) 

t (In.) 

Initial Condi 

tlons af Daslfa Varlablta 

(15.00) 

(310.00) 

(3500.0) 

(0.60) 

(3500.00) 

(0.40) 

1 

0.18491 

0.94687 

37.05 

445.34 

5834.00 

0.37 

6141.00 

0.13 

3 

0.30176 

0.098593 

39.81 

1019.00 

31000.00 

0.38 

5030.50 

0.33 

3 

1.0686 

1.3734 

30.93 

397.30 

1589.00 

0.43 

7837.50 

0.37 

4 

0.37654 

0.75959 

17.13 

435.00 

3655.50 

8.56 

4171.50 

0.31 

» 

0.43384 

1.1983 

31.67 

411.00 

1678.00 

8.64 

4136.50 

0.13 

6 

0.53838 

1.1863 

14.35 

403.30 

1696.50 

0.35 

4381.00 

0.38 


CASE 

08JECTIW 

FUNCTION(S) 

DESIGN 

VA«IA8tf(S) 

1 WalfNt 

t/N/l/Wt 

3 Cast 

1/H/l/Wt 

9 N 

1/H/l/Wt 

4 HilfNt/Cast 

8/Vt/VVt 

5 HalfUt/H 

VH/l/VVt 

6 Halftt/Cast/H 

1/H/L/VVt 
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Figure 1. Constrained to Unconstrained 
Conversion, Infeasable Initial Design x 



Figure 2. Constrained to Unconstrained 
Conversion, Feasable Initial Design x 



Figure 3. Constrained to Unconstrained 
Conversion, Constrained Minimum x Cm 



Figure 4. Constrained to Unconstrained 
Conversion with Multiple Objectives 
at Constrained Minimum 




Figure 5. Symetric Three Bar Truss 



0 20 40 to «o 
Figure 6. Weight/Cost Comparison, 
Three Bar Truss 



0 10 20 
Figure 7. Weight/Cost Comparison, 
Three Bar Truss 



Figure 8. Transmission Tower with Cables 



o 1 t 

Figure 9. Weight/Cost Comparison, 
Transmission Tower with Cables 
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